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1, Introduction 


In this semi-annual status report, the work previously - initiated 

dealing with the joint pdf which is a modification of the bivariate Gaussian 

pdf is discussed in detail and results are presented for a global reaction 

model using this joint pdf (Section 2). An alternative joint pdf, the ’’most- 

/2i 

likely" pdf, suggested by Pope ' ', is discussed in Section 3. A criterion 
is developed in Section 4 which permits the selection of temperature pdf's 
in different regions of turbulent, reacting flow fields. Section 5 sets forth 
two principal approaches to the determination of reaction rates in computer 
programs containing detailed chemical kinetics. These models are set forth 
in the context of current research in this field and, it is argued, represent 
a practical solution to the modeling of species reaction rates in turbulent, 
reacting flows. 



2, The Modified Bivariate Gaussian PDF 

Numerous experimental studies are available which describe pdf's for 
temperature in both non-reacting and reacting turbulent flows. A number 
of these (eg, "clipped" Gaussian, beta, ramp) have been discussed pre- 
viously. However, no experimental studies are available describing joint 
pdf's of temperature and species in any flow. The bivariate, or two di- 
mensional, Gaussian pdf has been chosen as a joint pdf which has at least 
a conceptual basis in physical reality. This conclusion is drawn by ima- 
gining the behavior of the marginal pdf's of temperature and species, each 
a one-dimensional Gaussian pdf, during the course of a reaction. The mar- 
ginal pdf's show that the mean temperature can be low while the mean con- 
centration of fuel, for example, is high. This could correspond to an 
initial stage of a reaction. As the reaction proceeds, the mean temperature 
increases as the mean concentration of fuel decreases. The marginal pdf's 
can exhibit this behavior. Based on this presumed physical behavior, it 
is assumed that the bivariate Gaussian pdf serves as a reasonable approxi- 
mation to the actual joint pdf of temperature and species in a turbulent 
reacting flow. This is intended as an Initial step in the formulation of 
such pdf's. 

The form of the bivariate Gaussian pdf is 
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The bivariate Gaussian pdf is valid on the interval of minus to plus in- 
finity for both variables and as such, cannot be utilized in its present 
form. To remecfy this situation, the pdf is modified as follows, Let 
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It follows that. 
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which is the form of the modified bivariate Gaussain pdf. This pdf is a valid 
one, as it satisfies the condition required of a pdf. 

Equation (6) below is a general expression for the reaction rate of 
fuel for any particular one-step, irreversible reaction of the form: 


F 0 O P 


from which 





where 
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for the fuel and the oxidizer. 
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That is, Xp and Xq are dimensionless concentrations that are bounded by zero 
and one. The terms Cp, Cq represent the initial concentrations of fuel and 
oxidizer, respectively,and n and m are constants. By assuming that n and m 
are both equal to unity, it is possible to rewrite (5) as 
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Therefore, Wp is a function of two variables, t and Xp. Equation (6) 
can be used along with the modified bivariate Gaussian pdf to determine the 


mean turbulent reaction rate. That is, 
1 1 
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where: 




Jj pct.x,) dt dXp 
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For a given or calculated mean turbulent reaction rate, the correspond- 
ing laminar term can be determined by inserting the proper values of t and 
x^, the mean dimensionless temperature and mean dimensionless concentration 
of fuel, respectively, into eq (6). The reaction rate amplification ratio 
is then: 




( 8 ) 
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This reaction rate amplification ratio considers the combined effects of 
temperature and species fluctuations on the reaction rate. As described in 
ref 1, Section 3(a), the assumption of statistical independence between tem- 
perature and species fluctuations leads to 



Amplification ratios calculated with both eqs (8) and (9) can be compared to 
determine the effect the inclusion of species fluctuations on the reaction 
rate, in addition to the temperature effect. For purposes of a parameter study, 
the following one-step, irreversible reaction is considered: 

^ i: *^2.^ 

The data of ref 3 are utilized to determine an approximate relationship between 
the dimensionless mean temperature and the dimensionless mean concentration of 
hydrogen in a turbulent diffusion flame. Shown in figure 1 are curve fits of 
data from ref 3, (fig 2, II 7 /U = 8). The mean temperature is non-dimension- 
all zed by use of the values of T_,^ = 2160K and 1 . - 300K. The range of 

iiiax nil n 

data for which the temperature is decreasing is ignored as it is assumed to 
be the result of heat transfer. In a hypothetical adiabatic system, this 
heat transfer would not occur. From the figure, it is readily determined 
that an approximate relationship between t and xj^ is: 

= 1 - t (10) 

This is deemed adequate for purposes of this parameter study. 
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Assume that in t*ie radial direction of a confined axi symmetric turbulent 


diffusion flame, the Prandtl mixing length hypothesis is valid. Thus; 



where L » mixing length and y ■ radial coordinate. 


Then 



1 at /dy 1 


1 


i’- 


The following data are extracted from ref 3, (Uj/Ug * 10, X/D * 40); 


r/a 



1-2 

0.08 

0.05 

2-3 

0.17 

0.07 

3-4 

0.24 

O.ll 

4-5 

0,22 

0.21 

5-6 

0.21 

0.23 

6-6.1 

0.30 

0.30 


Note: r/a y 
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Averaging yields 


* 0.16 

Then 



2 r-ir 

Thus, t' = 3/? has been chosen for the purposes of this parameter study. 

figure 1 can also be utilized for the determination of the correlation co- 
efficient, p, between temperature and species. The Statistics I Pac for the 
Hewlett Packard 41C programmable calculator was employed. This pac utilizes 
an approximate expression for thn correlation coeffecient. 

Pairs of data points for t and xfj^, at corresponding values of X/D up 
to approximately 130 were chosen. It was found that p - -0.99 is the value 
consistently calculated. It is assumed that this result is the same for in- 
stantaneous values of temperature and concentration. A similar procedure 
is follov/ad with the use of data from ref 3, and Identical results obtained. 

The results obtained using the modified bivariate Gaussian pdf, eqs (7) and '(8) 
are compared with the reaction rate amplification ratio of eq (9), determined 
by use of the beta pdf, in which only the effect of temperature fluctuations on 
the reaction rate is considered. This amplification ratio is equivalent to 
the reaction rate constant amplification ratio utilizing the beta pdf presented 
in refs I and 4. A comparison of these two reaction rate amplification ratios 
allows the combined effects of temperature and species fluctuations to be compared 
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with the temperature fluctuation-only effect, The cases considered are those 
which involve the variation of the following; mean dimensionless temperature 
and mean species concentration, dimensionless temperature and species concen- 
tration fluctuations, the correlation coefficient, the relationship between 

and and the maximum temperature. 

All numerical integrations were performed with the use of a Monte Carlo 
(random variable) technique. The number of random variables utilized is twenty- 
thousand (2 X 10*^). This number ensures that essentially no truncation error 
is incurred. There is a limit on the modified bivariate liaussian's range of 
applicability, generally in the ’region of low temperature and species concen- 
tration fluctuations, as well as high mean temperature and correspondingly low 
mean species concentration. This is due both to the mathematical behavior of 
the function and the formulation of the model. These limitations are reflected 
in the results which follow. 

Shown in Table 1. are the defining values and relationships for the cases 
considered. The initial values of concentration are taken from Test Case 1 of 
ref 5. All other constants are taken from ref 5 and are considered typical. 
Throughout this discussion, the use of the terms "temperature" and "concen- 
tration" denote mean dimensionless temperature and mean dimensionless species 
concentration, respectively. 

Shown in fig$ 2 through 6 are the results of the 



Table 1 Case Data 


Case 

No. 

B 

h 

(°K) 

T 

'max 

P 

vs. W 

1 

0 

10116 

2500 

-0.9 

XT’ « 2/3 
H2 

4 

0 

10116 

3000 

-0.9 

» 2/3 

”2 

5 

0 

10116 

2500 

-0.8 

^ - 2/3 

6 

0 

10116 

2500 

-0.9 

'% * 1/2 


Note; A * 8.4 X 10^®. 'f . « BOQOK, Cfi « 2.405 x 10'®gmol/cni\ 

iMin ng 

C §2 * 1*5 X 10“^gn]ol/cm% » 1 - t, for all cases 

variation of the reaction rate amplification ratio with combined temperature 
and concentration fluctuations, at various constant values of temperature and 
concentration. Selected results of the temperature-only case (beta pdf) are 
shown for purposes of comparison, 

Figs 2 and 3 show that at high temperatures and low concentrations, the 
amplification ratio of the combined case is greater than that of the temperature- 
only case. Fig 4 indicates that at moderate values of temperature and concentra- 
tion, there are values of the fluctuations that cause the amplification ratio 
of the combined case to become less than that of the temperature-only case. 

This trend is continued in fig 5, where the amplification ratio of the former 
15 always less than that of the latter. In all these figures, both cases show 
an increasing amplification ratio with increasing fluctuations. In fig 6 for 
t = 0.6, the amplification ratio of the combined case at first increases, reaches 
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a relative maximum, then decreases as the fluctuations increase . With 
t * 0,7, the amplification ratio of the combined case is less than unity 
and decreases with increasing fluctuations. The amplification ratio of the 
temperature-only case is greater than that of the combined case for both values 
of T and increases with increasing fluctuations. Inspection of all these fi- 
gures shows that the amplification ratios of both cases decrease with increa- 
sing temperature and a corresponding decrease in concentration, at constant 
values of fluctuations. The following is a postulated explanation for the 
behavior of the combined case amplification ratio. 

As previously shown (refs 1 and 4), the effect of temperature fluctuations 
is to amplify the Arrhenius reaction rate constant. The consideration of the 
effects of concentration fluctuations on the reaction rate involves two fac- 
tors which are in addition to this temperature effect. These factors are the 
probability of collision between molecules of fuel and oxidizer and the time 
available for two adjacent molecules of these types to react, A proper com- 
bination of, these two factors allows thi reaction to occur. 

At low values of temperature and correspondingly high values of concentra- 
tion, the fluctuations of concentration provide an amplification effect that 
is in addition to the amplification effect of the temperature fluctuations. This 
effect increases as the fluctuations increase (fig 2), This may be interpreted 
in terms of a higher effective concentration. While the molecules of fuel and 
oxidizer are spending less time, on the average, in individual collision vol- 
umes, there is a high enough concentration of them so that the probability of 
collision is very high, on the average, and leads to a resultant high overall 
probability of reaction. Thus, as a result of the increasing concentration 
fluctuations, the probability of collision increases more than the probability 
of molecules spending enough time in their collision volumes decreases. As 
the temperature increases and the concentration decreases, this effect is re- 
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duced (fig 3). That is, the probability of collision is not as large as at 
higher concentrations. Up to this point, the amplification effect of the con- 
centration fluctuations is still in addition to that of the temperature fluctu- 
ations. Hence, the amplification ratio of the combined case is greater than 
that of the temperature-only case. The amplification ratios of both cases in- 
crease with increasing fluctuations. At some point, the value of concentration 
is not high enough, and an increase in the fluctuations causes the amplification 
ratio of the combined case to become less than that, of the temperature-only case 
(fig 4). This is because the effect of the concentration fluctuations at this 
value of concentration is to cause a very low probability of two molecules spend- 
ing enough time in the required collision volumes. This is not sufficiently 
countered by the increased probability of collision caused by the concentration 
fluctuations and temperature amplification effect, However, the temperature am- 
plification effect is still strong enough to cause the amplification ratio to 
increase with increasing fluctuations. This is further demonstrated in fig 5. 
Eventually, the concentration and probability of collision become small enough 
so that at some value of the fluctuations, the amplification ratio decreases 
with increasing fluctuations (fig 6). Thus, at this point, the effect of the 
temperature fluctuations is not large enough to cause an increase in the ampli- 
fication ratio with increasing fluctuations. Finally, the concentration is so 

I 

small, the effect of the concentration fluctuations is to cause the mean turbu- 
lent reaction rate to become less than the corresponding laminar value. Hence, 
an amplification ratio less than unity results (fig 6). Therefore, in a given 
region of a flow, there are so few molecules of fuel and oxidizer, the probabi- 
lity of any two of them spending enough time in their collision volumes is very 
small, as is the probability of collision. This combines, despite the tempera- 
ture amplification effect, to give an overall probability of reaction which is 
lower than the corresponding laminar value. 
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The situation of a mean turbulent reaction rate which is less than the 
corresponding laminar valus has been observed through a combination of analysis 
and experiment in a turbulent premixed flame. This is reported in ref 6. In 
that case, calculation of the mean turbulent reaction rate through the solution 
of the fuel species probability transport equation results in much better agree- 
ment with experiment than the calculated average temperature laminar reaction 
rate. The actual mean turbulent reaction rate Is less than the corresponding 
laminar value and the use of the probability transport equation predicts this 
to be so. 

Shown in fig 7 is the variation of the amplification tatio with the maxi- 
mum temperature, for the case of combined temperature and species fluctuations. 

As seen from the figure, a change in the maximum temperature can cause a very 
substantial change in the amplification ratio. An increased inaximum tempera- 
ture causes a decrease in the amplification ratio. This behavior is similar 
at various values of t and xi. '. 

H2 

In fig 8 is the variation of the amplification ratio with the correlation 
coefficient, for the case of combined temperature and species fluctuations. As 
shown, an approximate increase of twelve percent (-0.9 to -0.8) causes only a 
very small change in the amplification ratio. The change which does occur is a 
decrease in the amplification ratio. This behavior is similar at various values 
of t and x^. Figs 9 and 10, comparing p = -0.9 and p = -0.5, further illustrate 
the relative lack of sensitivity of Z' to the correlation coefficient. 

Fig 11 shows the variation in the amplificeition ratio with a change in the 
relationship between and for the case of combined temperature and species 

fluctuations. As demonstrated, an approximate change of thirty percent (2/3 to 1/2) 
causes only a slight change in the amplification ratio. In general, the reduced 
concentration fluctuations result in a smaller amplification ratio. This be- 


ll 
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havior i3 similar at various values of t and 

The results represented by Figs 2 through 11 may be summarized by the 
following conclusions: 

1. The reaction rate amplification ratio of the combined case 
may be greater than, equal to or less than, that of the 
temperature-only case. 

2. The amplification ratio of the combined case increases with 
increasing temperature and species concentration fluctuations, 
at low values of temperature and correspondingly high values 
of concentration. At high values of temperature and corre- 
spondingly low values of concentration, it decreases with 
increasing fluctuations. Thus, the model employing the mo- 
dified bivariate Gaussian pdf predicts that at an initial 
stage of a reaction (low temperature and high concentration) 
the combined case mean turbulent reaction rate will be great- 
er than the corresponding laminar value. As a reaction pro- 
ceeds, this mean turbulent reaction rate will become less 
than the corresponding laminar value. 

3. The amplification ratio of the combined case shows some sen- 
sivity to changes in the value of maximum temperature. Hence, 
this value must be carefully selected. 

4. The amplification ratio of the combined case is relatively in- 
sensitive to : (a) small changes in the correlation coefficient, 

(b) moderate changes in the relationship between temperature and 

s 

species fluctuations. 


12 



3. The "Most-L1kely" PDF 


A joint pdf alternative to the modified bivariate Gaussian discussed in 
Section 2 is that suggested by Pope and is termed the "most-likely" pdf. 
The single-variable form of this pdf is 

U) ’ + X ^ t’’ j ( 11 ) 


where the X*s are obtained from the following constraint equations: 

fpCtlJt.l 

o 

t * f ‘ p (13) 

o 

1 

' ((t-t)^p^t)dt (14) 

0 

Hence, knowledge of, say, the mean dimensionless temperature and the mean square 
temperature fluctuations permits a determination of X^, X, and X 2 using eqs (12), 
(13) and (14). Currently, this one-variable form is under investigation for com- 
parison with the results obtained using the beta and ramp pdf's. 

Multi-variable forms of the most-likely pdf are: 


Two-variable 


p(x, [x. <- X,x, f X,x,x^ j 


(15) 


Three-variable 


f X^ %,X3 + X4 ] 


(16) 
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In eq '15), thru X 3 are determined from the constraint equations; 

1 1 


1 1 
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X 


1 



(5(x,,xjdx,<ix^ 


%l X 


/ 

X . 


1 1 



X, )(x-X^) . 
p(x,,x,)dx,</Xj^ 


(17) 

(18) 

(19) 

( 20 ) 


In eq (16), X thru Xs are determined from equations equivalent to eqs (17) 
thru (20) for~,T-, x . xU' , x'x'. x'x', and 

1 2 3 i 2 ‘ ] 3 2 3 

1 < 1 

( 5 \ 3 -Xj)" 1 

o o o 


The term (j^in eqs (15) and (16) is an "a priori" probability whose two- 
variable form is 



( 21 ) 
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where X is the ratio of the time scale of the turbulence to the time scale 
of the reaction; ie, 


X 



( 22 ) 


where * K]e, where k: is the turbulence kinetic energy and c its dissipation 
rate. A determination of the reaction time scale, is less clear-cut. In 
ref 20, Pope's selection is arbitrary and unexplained. One way of formulating 
this term at a particular point in a turbulent flow field might be to construct 
an "effeccive half-life time" at that point from 



(23) 


where k is the rate constant obtained by averaging the rate constants previously 
obtained at neighboring points in the flow and "c is an average molar concentra- 
tion obtained at those neighboring points. For example, for the bimolecular 
elementary reaction A + B = C + D, c = (c^ where c^ and are turbulent 

mean values obtained at neighboring points. 

The covariance, eq (20), may be obtained either from a transport equation or, 
more simply, from the relationship 


x'x' 

1 i 


= p(x^)^ (^^) 


(24) 


Where p is, as m Section 2, the correlation coefficient. Since the results ob- 
tained using the modified bivariate Gaussian were relatively insensitive to the 
value of p, a value of p = 0.^. is recommended as an initial estimate. 
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The formulation discussed above requires, therefore, a knowledge of first 
and second moments; that is, x x^ and along with an estimate of the 

parameter X in eq (21). The manner in which this analysis would be incorporated 
into a large-scale computer program describing turbulent reacting flows is dis- 
cussed in Section 5. The three-variable fonnulation closely parallels the two- 
variable problem outlined above and will be discussed in detail in the next Status 
Report. 

4. A Criterion for Selecting a Temperature PDF in Turbulent Flows 

On the basis of comparitive studies which represented turbulent temperature 
fluctuations by "clipped" Gaussian, beta- and temperature-ramp pdf's (refs. 1 and 
4), it was concluded that amplification ratios are substantially higher in regions 
exhibiting ramp-like behavior. 

It is anticipated that in complex turbulent flow fields, regions will frequent- 
ly exist in which rainp-like temperature fluctuations will be present. Outside 
these regions, temperature fluctuations may be best represented by, for example, 
the beta pdf. Hence, a criterion is required for selecting an appropriate tempera- 
ture pdf when making calculations describing turbulent reacting flows. Such a 
criterion is developed in this Section. 

Temperature fluctuations which exhibit a ramp-like structure have been ob- 
served experimentally under various flow conditions. Fiedler studied a 
heated two-dimensional mixing layer, i.e,, the higher velocity stream was heated 
while the external stream had zero velocity and was at ambient temperature. A 
turbulent heated jet with a coaxial flowing' stream was investigated- by Antonia, 
Prabbu, and Stephenson. La Rue and Libby and Gibson, Chen and Lin 
independently studied the turbulent flow behind a heated cylinder. Ramp-like 
temperature fluctuations were observed in all four investigations. Despite the 
fact that all these experiments were performed in non-reacting flows, it is assum- 
ed that these results can be generalized to reacting flows since ramp-like fluct- 
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uations were observed under a wide variety of conditions likely to be encountered 
in reacting flows. Gibson et al. believed these types of fluctuations are due 
to sharp thermal gradients, with the point of highest temperature of the ramp 
occuring at either the downstream or upstream end of the turbulent/irrctational 
interface depending on the sign of the vorticity of the main flow. 

Experimentally obtained temperature signals are shown in fig 12 
Note the similarity of the ramp-like structure of the signals even though they , 
were obtained under different flow situations, 

Fig 13 shows data obtained by Fiedler The ordinate is the flatness 

factor F (also called Kurtosis) and the abscissa is a normalized distance. Note 
that there is a horizontal line through F=3 labeled "Gaussian distribution," 
since the Gaussian pdf has the singular value of 3 for the flatness factor. It 
possesses a singular value since there are no adjustable constants, as there are 
in the ramp pdf. 

The general trend of the flatness factor in fig 13, is that it has a value 
of approximately 3 at the centerline (h = 0), decreases to a value near 2, then 
sharply increases in the outer edge of the flow field. The same trend was ob- 
tained by Antonia et al. as shown in fig 14. The distance parameter is not 
the same as Fiedler's due to the use of different normalizing constants. Note 
the similarity between the figs 13 and 14. This trend of the flatness factor 
was observed at different axial distances and for different velocity ratios. 

This general behavior led Fiedler to hypothesize that the value of 2 was charac- 
teristic of the "sawtooth" appearance of the temperature signals . His criterion 
was supported by the appearance of ramp-like temperature signals at values of 
n near the centerline. For positive values of n (n = 0,0.5), the signals exhibit 
a definite ramp-like structure. For these values of n» the value of the flatness 
factor lies between 2 and 3. Antonia et al. observed the same phenomenon. 
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Clearly then the flatness factor F can serve as the criterion for specify- 
ing when the temperature-ramp pdf should be used in a given region of a tur- 
bulent flow field, In order to specify a range of F (around F * 2.0) in which 
ramp-like temperature fluctuations will be present, an analysis will be developed 
to determine the moments of the ramp pdf. 

It will be recalled, refs 1 and 4, that the temperature-ramp pdf is 

«• 

U*)- f ' (25) 

where an asterisked quantity is one which has been normalized by dividing by a 
constant A*, an upper-case quantity is a dimensionless, temperature-related 
function, and a lower-case quantity is a particular value of the function. The 
ramp portion of the fluctuations is represented by Pj^.^ (r*): 



and the Gaussian portion isj 




-ti^r 


where c, g and are constants. The moments are given by 

•e 

-00 ' 

Substituting eqs (25), (26) and (27) into (28) leads to 

E (t*'') * f ft*-'’ ^ 

' 1 } ./liT <r* ^ . 
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which is valid for n « 1, 2, 3 and 4. The gainma function In eq (30) can be re- 
placed by 




As a result of eqs (30) and (31), eq (29) now becomes 



where the integration limits have been changed to reflect the range of r*. 

A Monte Carlo random variable integration technique was employed to carry 
out the required integrations. Results for the flatness factor, F, are shown 
in figs 15 through 18 for the following ranges of the constants: 

1,1 i k‘ < 1.^ 

0.1 4 <r*4 0.5 
0.75 4 t i I.IS 
1,0 4 A 4 1.0 
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The values recommended by Antonia and Atkinson are «* * 1.8, a* ■ 0,25, 
c » 1,0 and 6 » 1,25. Arbitrarily selecting values around these of ±5%, the 
recommended range for the flatness factor Is: 

2.0 < F < 2.4 

wmt mmm 

That is, within this range of F, temperature fluctuations should be represented 
by the temperature- ramp pdf. Outside this range of F, the beta pdf or some 
other appropriate pdf should be used to represent the temperature fluctuations, 

The manner In which this criterion might be incorporated Into a large-scale 
computer program testing reacting turbulent flows is detailed In Section 5 herein. 

5. Discussion of Alternative Models for the Determination of Reaction Rates in 
Turbulent Reacting Flows 

In order to place into context the models proposed in this Section, the five 
principal approaches of current interest will initially be discussed. The approach 
es referred to deal with the determination of species reaction rates In a turbu- 
lent reacting flow. 

The first approach involves an expansion of the governing conservation equa- 
tions of overall mass, momentum, species and energy In terms of limiting values 
of certain parameters. This is described in ref 12. The status of this method 
has advanced rapidly in recent years, but the resultant expansions are often found 
to be perturbations of laminar flows and thus may only apply to limited cases of 
turbulent flows. 

The second approach simply utilizes the mean temperature (average, static or 
laminar) for the calculation of the Arrhenius reaction rate constant of the parti- 
cular reaction under consideration. This is coupled with a modeled effect of the 
"unmixedness" on the reaction rate. The term "unmixedness" refers to the situa- 
tion In which the turbulent fluctuations of fuel and op’dizer species are so great 
that reaction ceases at various points in a flow. This determination of a turbu- 
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lent reaction rate is utilized with existing turbulence modeling and is describ- 
ed in ref 5. While not entirely neglecting the effects of the species fluctua- 
tions. this approach does neglect the effect of the temperature fluctuations and 
other possible effects of the species fluctuations on the reaction rate. It has 
been shown in ref 13 by means of kinetic theory that the effective activation 
energy of any reaction in a turbulent flow is lower than that of the same reac- 
tion in an equivalent laminar flow (same mean temperature). This implies that 
the effect of the turbulent temperature fluctuations is to present to the flow 
a higher effective temperature and amplified Arrhenius reaction rate constant. 

The third approach models the turbulent reaction rate based on observed 
physical phenomena in a flow. The result of the application of this method is 
termed the "eddy break up" model and is described in ref 14. This method lacks 
rigor and may be limited by an incomplete understanding of the physical pheno- 
mena Involved. 

* 

The fourth approach utilizes the concept of the probability of temperature 
and species as being transportable quantities in a flow. The transport processes 
are those of molecular diffusion, turbulent mixing, etc. Thus, a transport equa- 
tion, utilizing existing turbulence modeling, can be written in a general form 
and solved for specific flows. The resulting probability distribution can then 
be used in conjunction with the tools of applied probability theory to calculate 
the properties of a turbulent reacting flow. The modeling process is described 
in ref 15 and a transport equation solution procedure is described for example 
in ref 16. This method is attractive from a fundamental point of view but un- 
fortunately, can require large amounts of computer time and storage. This pro- 
blem arises when it is necessary to consider multi step reaction mechanisms, as 

s 

is often the case. This is seen in refs 17 and 18, in which predictions of 
nitric oxide formation in the turbulent exhaust from gas turbine combustors are 
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made and in ref 5» where it Is ne^'essary to consider a twenty-five elementary 
reaction mechanism for the supersonic combustion of hydrogen with oxygen. 

The fifth approach is similar in principle to that just described and 
has been the subject of the present study. However* instead of solving a 
transport equation for probability, a specific probability density function 
(pdf) is assumed in particular regions of a turbulent, reacting flow field and 
the mean turbulent reaction rate is then determined. This method is attractive 
as a solution to the probletn of turbulent reacting flow prediction because of its 
basis in a fundamental physical concept and its relative computational simplicity. 

Within the framework of this "presumed pdf method" three models (at 

least) may be formulated for the reaction rate; 

I. A model which assumes statistical independence between temperature and 
species fluctuations, in which temperature fluctuations are represented by appro- 
priate pdf's in the various regions of the flow and the effects of species fluct- 
uations, are treated using "unmixedness" criteria or some other equivalent 
formulation. 

II. A three-variable model which assumes statistical independence between 
temperature and species fluctuations, but in which the latter are treated using 
a joint pdf. Thus, for example, for a reaction rate term of the form; 


this model yields a mean turbulent reaction rate given by; 

. j KCt) pCt^t j . 


( 33 ) 


where it has been assumed that c, 


0 , min 
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Ill, A three-variable model using a joint pdf for the temperature and 
species concentrations. The form of the reaction rate tews In this case Is; 


W. a C,, 

1 


1 1 1 
®®0 


The latter two models mentioned above will now be discussed In detail. 

In model II, the effect of the temperature fluctuations on the rate constant 
is treated using the methods discussed herein and In previous Status Reports 
\ In most regions of the flow, the beta pdf, for example, would be used. In 
regions wherein the flatness factor Is in the range from 2.0 to 2.4, the ramp 
pdf should be used. By analogy with the temperature fluctuations, 
the flatness factor may be determined from a transport equation of the form 

II: i 1/ \..i 


aix. 


V 






sy - L a) J - ^ 

where 0 = t/(T^)\ Assuming local equilibrium exists, the generation and dissi- 
pation terms on the right-hand-side may be set equal, leading to 




where the equality p = ^ has been used. It will be recalled that 

can be obtained from a similar relationship: 


C. W it 

> VTf j 
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The constants c i c , c and c* must be determined empirically. 

12 1 4 


The joint pdf for the species In eq (33), p (x^^, x^^), might be taken to 
be of the modified bivariate Gaussian form (Section E) or of the "most-likely” 
pdf form (Section 3). That 1s, 


‘•‘P (38) 

where cj , x-, Xi, Xg, and Xj are obtained using the methods described In Section 
1 - Jierel.i. Clearly, for a chemical system of some complexity, these methods would 
require Inordinate quantities of computer time. A preferable procedure would be 
to obtain in advance a tabular array or polynomial curve fits from which value of 
the coefficients of p can be obtained by inputing x^^, 'x^^, *i^, xj^, and X (the 
ratio of turbulent- to-reaction time scales). Then, using data from adjacent n'odes, 
estimates of the input parameters can be used to enter the table or to use the 
curve fits to obtain ^2 ^ 3 * The reaction rate terms are next cal- 

culated using eq (33) for each term*, then the reaction rate Itself Is determined 
from the summation of the Individual terms. That Is, for a set of bimolecular 
elementary reactions of the general form 

the laminar reaction rate of, say, species Aj is 


% ‘ ■^“=A,‘=A. ^•-.<^A.'A4--- 


and the turbulent reaction rate is 
'^Ai* "^‘*Aj "••• 
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•ym M M 

Where w-,* is given by eq (33) and the remaining terms, w • /. , w a , etc. are 

1 nj "1 2^1 

given by expressionsequivalent to eq (33). Here again, the concept of using a 
tabular array or curve fits will serve to eliminate the necessity for repeated 
integrations of eq (33). That is, after selection of the appropriate p (t) using 


the flatness factor criterion, a table will be entered with T, )T^^, 

xP' and X to obtain 2-'; - ^ ^ 

”2 ^ ^ -I r / / 


7 




[(K,Ct)pCtWt 


Since w a * -k i. is readily determined using existing procedures w * 

1,^1 I ni M2 l»^l 

and like terms can be calculated. 

A major near-term task under this grant is the determination of values of 
■2-Mt, t^, x^^. x^^, )^, X). 

The third model mentioned earlier is a three-variable model in which a joint 
pdf is employed for both temperature and species fluctuations; ie, p(t, x^^, x^^^j 
In this case, the pdf will be assumed to be given by the following three-variable 
form of the "most-likely pdf": 

'Xjj.t Xa, X(.'>-A,x^^j (40) 


where the method for obtaining q and the constraint equations for X^, Xi, ...» Xe 
are similar to those detailed in Section 3 herein. As with model II above, the 
amplification ratio 2-' « = w^ «-/Wi « can be determined a priori as a function of 
T, 7^^, 7^^, t'-x^, t ' x^^ , x^x^^ and X and set up in the form of a tabular array or 
curve fits for use in large-scale computer programs. 
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It will be noted that the models discussed above require values of the mean 
square species fluctuations, , as well as the mean square temperature fluctua- 
tions, As is the care with the latter, transport equations for the species fluc- 
tuations could, conceptually, be introduced for each species as was done in ref 5 
for the fuel and oxidizer. The general form is 



where the symbols are the same as in ref 5 and In actual practice, 

computer time and space restrictions will generally militate against the addition 
of these differential equations, particularly when detailed chemical kinetics are 
involved. In light of the other uncertainties in the analysis, it is argued that 
a local equilibruim assumption is in order, similar to that used to obtain eq (36), 
in which case eq (41) reduces to 



Attention will be focussed during the forthcoming semi-annual period on the 
Model II problem discussed above. In addition, a one-dimensional, temperature- 
only form of the "most-likely" pdf^: 

will be investigated for comparison with the beta and ramp pdf models. 
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LIST OF SYMBOLS 


A 

a 

8 

C 

c 

F 

f (x,y) 
g 

k (t) 

k k* 

1 2 

L 

P (t) 

P (x,y) 

q 

r 

T 

t 

U 

U, V 
V 

w 

X 


pre-exponential factor, eq (7) 
constant used in eq (25) 
duct radius 

temperature exponent, eq (7) 

molar concentration; constant 

constant; molar concentration; constant in eq (26) 

flatness factor or Kurtosis 

two-variable pdf 

rms species fluctuation 

Arrehenius reaction rate constant 

constants 

Prandtl mixing length 
Prandtl mixing length 
one-variable pdf 
two-variable pdf 
"a priori" probability, eq (21) 
radial coordinate 
temperature 

"activation temperature" ratio of activation energy 
to gas constant 

dimensionless temperature 

axial velocity 

two-dimensional or axisymmetric velocities 

defined in eq (2) 

reaction rate i 

ratio of turbulence-to-reaction time scales, eq (22 
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I 


X 


X, y 
1 

«' 

3 

e 

X 

K 

P 


dimensionless molar concentration 
two-dimensional or axisymmetric space coordinates 
rate constant "amplification ratio," eq (9) 
reaction rate "amplification ratio," eq (8) 
constant in eq (26) 
constant in eq (26) 

dissipation rate of the turbulence kinetic energy 
parameters in the "most-likely" pdf 
turbulence kinetic energy 

« 

correlation coefficient; mass density 


u 

0 

0 

T 


Gaussian mean 
Gaussian rms fluctuation 
defined after eq (35) 
time 


Subscripts 

A , 

■ 

F 

I 

max 

min 

0 

R 

T 

t 

X 


species index 
fuel 

laminar value 
maximum value 
minimum value 
oxidizer 
reaction 
turbulent 

dimensionless temperature; turbulent value 
dimensionless molar concentration 
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Superscripts 

0 initial value 

* quantity normalized by dividing by A^. 
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Figure 2.. Variation of the reaction rate amplification 

ratio with temperature and species fluctuations 
for the combined (modified bivariate Gaussian 
pdf) and temperature-only (beta pdf) cases 
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Figure 3 . Variation of the reaction rate amplification 

ratio with temperature and species fluctuations 
for the combined (modified bivariate Gaussian 
pdf) and temperature-only (beta pdf) cases 
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Figure 4 . Variation of the reaction rate amplification 

ratio with temperature and species fluctuations 
for the combined (modified bivariate Gaussian 
pdf) and temperature-only (beta i)df) cases 
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Figure S . Variation of the reaction rate amplification 

ratio with temperature and species fluctuations 
for the combined (modified bivariate Gaussian 
pdf) and temperature-only (beta pdf) cases 
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Figure (# t Variation of the reaction rate amplification 

ratio with temperature and species fluctuations 
for the combined (modified bivariate Gaussian 
pdf) and temperature-only (beta pdf) cases 
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Figure 7 . Variation of the reaction rate amplification 
ratio with the maximum temperature for the 
combined case (modified bivariate Gaussian pdf) 
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Figure d . Variation of the realction rate amplification 

ratio with the correlation coefficient for the 
combined case (modified bivariate Gaussian pdf) 
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Variation of the reaction rate amplification 
ratio with the relationship between temperature 
and species fluctuations for the combined case 
(modified bivariate Gaussian pdf) 
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